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The production/absorption rate of particles in compressed and 
heated asymmetric matter is studied using a Relativistic Mean 
Field (RMF) transport model with an isospin dependent collision 
term. Just from energy conservation in the elementary produc- 
tion/absorption processes we expect to see a strong dependence 
of the yields on the basic Lorentz structure of the isovector ef- 
fective interaction, due to isospin effects on the scalar and vector 
self-energies of the hadrons. This will be particularly evident for 
the ratio of the rates of particles produced with different charges: 
results are shown for it + /-k~, K + /K° yields. 

In order to simplify the analysis we perform RMF cascade sim- 
ulations in a box with periodic boundary conditions. In this way 
we can better pin down all such fine relativistic effects in particle 
production, that could likely show up even in realistic heavy ion 
collisions. In fact the box properties are tuned in order to reproduce 
the heated dense matter formed during a nucleus-nucleus collision 
in the few AGeV beam energy region. 

In particular, K +, ° production is expected to be directly related to 
the high density behaviour of the symmetry energy, since kaons are 
produced very early during the high density stage of the collision 
and their mean free path is rather large. We show that the K + / K° 
ratio reflects important isospin contributions on the production 
rates just because of the large sensitivity around the threshold. 
The results are very promising for the possibility of a direct link 
between particle production data in exotic Heavy Ion Collisions 
(HIC) and the isospin dependent part of the Equation of State 
(EoS) at high baryon densities. 
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1 Introduction 



Recently the development of new heavy ion facilities (radioactive beams) has 
driven the interest on the dynamical behaviour of asymmetric matter, see 
the recent reviews [1,2]. Indeed, the isospin-dependent part of the EoS is of 
crucial importance in extrapolating structure calculations beyond the valley 
of stability and in astrophysical processes, such as structure and cooling of 
neutron stars and supernovae explosions. 

Here we focus our attention on relativistic heavy ion collisions, that provide a 
unique terrestrial opportunity to probe the in-medium nuclear interaction in 
high density and high momentum regions. We will show that particle produc- 
tion/absorption (here pions and kaons) processes in a dense and hot neutron 
rich medium are nicely sensitive to the details of the relativistic structure of 
the effective interaction. The results appear very promising for the possibility 
of directly pin down from data such fundamental microscopic information. 

Within a covariant picture of the nuclear mean field, for the study of asym- 
metric nuclear matter certainly the same considerations hold for the density 
behaviour of the symmetry term as in the iso-scalar case. The equilibrium con- 
ditions are directly related to the Lorenz structure of the iso-vector effective 
fields in a similar way as in the iso-scalar case (competition effects between 
attractive scalar and repulsive vector fields). However, for the description of 
the a± parameter of the Weizsaecker mass formula (in a sense equivalent to 
the ai parameter for the iso-scalar part), extracted in the range from 28 to 36 
MeV, there are different possibilities: (a) considering only the Lorentz vector 
p mesonic field, and (b) both, the Lorentz vector p (repulsive) and Lorentz 
scalar S (attractive) effective fields [3,4]. The second assumption corresponds 
to the two strong effective a (attractive) and uj (repulsive) mesonic fields of 
the iso-scalar sector. In both cases one can fix the empirical value of the 
parameter, however important effects at supra-normal densities appear due to 
the introduction of the effective 8 field. In fact the presence of an isovector 
scalar field is increasing the repulsive p-meson contribution at high baryon 
densities [4,2]. This is a pure relativistic effect, due to the different Lorentz 
properties of these fields (the vector p field grows with baryon density whereas 
the scalar 5 field is supressed by the scalar density). As a consequence, the 
isospin-dependent part of the EoS becomes "stiffer" when both fields are ac- 
counted for. This important feature appears in several models, see ref. [5] 
where the covariant structure of the symmetry energy was investigated adopt- 
ing different approaches to asymmetric nuclear matter such as Non-Linear 
RME, Density Dependent Hadronic (DDH) and Dirac-Brueckner-Hartree- 
Fock (DBHE) theories. Moreover, relativistic approaches to asymmetric nu- 
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clear matter naturally lead to an effective (Dirac) mass splitting between 
protons and neutrons due to the Lorentz scalar nature of the 8 field, [3,4,2]. 
This is of relevance for the dynamics of heavy ion collisions and it has recently 
been microscopically studied within the DBHF approach [6,7]. 

The study of the influence of the high density symmetry energy on the HIC 
dynamics has been recently started. So far one has considered collective flow 
observables, such as in-plane and elliptic flows [8-11] for protons and neutrons 
(or light isobars), isospin equilibration (or transparency) and pion production 
in heavy ion collisions at intermediate relativistic energies up to 1 — 2 AGeV 
[12,5,13]. However, definite conclusions have not been drawn yet due to either 
the lack of experimental information on some observables or the moderate 
dependence on the density behaviour of the symmetry energy. Comparisons 
with recent experimental data on isospin transparency seem to suggest a stiff 
high density behaviour of the symmetry energy, as predicted by RMF theory 
with the inclusion of the 5 meson [13,14]. 

In order to explore the symmetry energy at supra-normal densities one has to 
select signals directly emitted from the earlier non-equilibrium high density 
stage of the heavy ion collision, possibly without any secondary interaction 
with the hadronic environment. 

In fact particle production is of particular interest since it mainly takes place 
in the first binary collisions during the formation of the high density phase, 
where the energy densities are still above resonance thresholds. This has stimu- 
lated recent experimental studies [15,16]. However, the investigation of particle 
production represents a non-trivial task due to the following reasons: 

- (a) At intermediate energies (< 1.6 AGeV) some of the pions created during 
the high density phase are directly emitted, but many of them undergo 
secondary re-scattering with isospin exchange before emission, [5]. 

- (b) Production and reabsorption processes must be consistently implemented 
in relativistic transport models with isospin effects. In fact within an isospin 
dependent mean field, isospin exchange processes such as R' <-> nN, irirR 
(N,R,R' stand for protons/neutrons and the A ±,0 ' ++ isospin resonance 
states, respectively, and 7r ± ' for pions) modify the scalar and vector con- 
tent of the available energy between the in- and outgoing channels. This has 
to be taken into account in the energy conservation, and it is particularly 
relevant around the thresholds. 

- (c) The situation for the production of particles with strangeness such as 
kaons (K ±,Q ) is similar, but even more complicated due to the availability 
of more than 50 isospin channels for the creation of K + and K°, only from 
BB and Bn processes (B stands for a resonance or a nucleon). 

Kaons are produced very early during the high density stage of the collision, 
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see the compressibility dependence in refs. [17,18]. Moreover, the production 
of K +, ° is expected to be directly related to the high density behaviour of the 
symmetry energy, since the mean free path of these low energy kaons in nu- 
clear matter is rather large (X(K + '°) ~ 7 fm [17,18]). These two facts make the 
subthreshold kaon production a very promising probe for the high density sym- 
metry energy. Anti-kaons (K~) do not fulfill the desired criteria due to their 
relatively small mean free path value (strong reabsorption effects). A related 
problem for K~'s is the importance of off-shell contributions in the propaga- 
tion and the corresponding uncertainty on the self-energy evaluations [19]. At 
variance the quasiparticle approach followed here in a transport RMF frame 
is consistently justified for K + '° mesons. In any case the anti-kaon production 
is not so relevant at the energies of interest in the present study [20,21]. 



In this article we investigate hadronic matter properties in a finite box, with 
periodic boundary conditions, performing cascade calculations by an extended 
isospin dependent collision integral. The temperature and baryon density con- 
ditions are tuned in order to reproduce the early stage of a realistic heavy ion 
collision at beam energies around the threshold for kaon production. Varying 
the asymmetry of the system we can then study the influence of the isovector 
part of the hadronic in medium interaction. 



Thus, we discuss particle production in fixed conditions of density and temper- 
ature (for both symmetric and asymmetric matter). In this way we can have 
a precise control on the production/absorption rates in the different channels. 



The description of the mean field is important, since nucleons and resonances 
are dressed by the mean field self-energies. This will directly affect the thresh- 
old energies, in particular for inelastic channels. We first give an outline of 
the mean field models of asymmetric nuclear matter, before passing to the de- 
scription of the box calculations (Sect. 2). The modification of the threshold 
conditions according to energy-momentum conservation is a non-trivial task 
which will be discussed in Sect. 3. The isospin dependence of the chemical 
equilibrium conditions, to which particle multiplicities obey after a given time 
interval t eg , is illustrated in Sect. 4, considering in particular the n~ /n + ratio. 
The results are presented in Sect. 5. We will first discuss isospin effects on 
resonances and pion production, before passing to kaons. The results provide 
an important guide to the understanding of particle production in heavy ion 
collisions at relativistic energies. 
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U (fm 2 ) 


U ifm 2 ) 


f P (fm 2 ) 


fs (fm 2 ) 


A (fm- 1 ) 


B 




NL 


9.3 


3.6 


0.0 


0.0 


0.015 


-0.004 




NLp 


9.3 


3.6 


1.22 


0.0 


0.015 


-0.004 




NLp5 


9.3 


3.6 


3.4 


2.4 


0.015 


-0.004 





Table 1 

Coupling parameters in terms of fi = (^-) 2 for i = a, u, fc = (^-) 2 for % = p, 5, 
A = and B = £ for the non-linear NL models [23] using the p (NLp) and both, 
the p and 5 mesons (NLp5) for the description of the isovector mean field. The NL 
model does not contain any isospin dependence. 

2 Cascade calculations for hadronic asymmetric matter 

2.1 Relativistic effects on the symmetry energy 

Within a relativistic framework the energy conservation in binary collisions 
between hadrons in nuclear matter is directly related to the Lorentz com- 
ponents of the hadron self-energy E = E s — 7 M S M . The scalar and vector 
self-energies E s , E M generally depend on baryon density and momentum, ac- 
cording to microscopic DBHF calculations [22]. However, since we are focus- 
ing here on isospin effects, i.e. on the iso-vector part of the self-energies, and 
DBHF studies of asymmetric nuclear matter are still rare or just being started 
[6,7], we will use a simple phenomenological version of the Non-Linear (with 
respect to the iso-scalar, Lorentz scalar a field [23]) Walecka model which 
corresponds to the Hartree or Relativistic Mean Field (RMF) approxima- 
tion within the Quantum-Hadro-Dynamics [24]. According to this model the 
baryons (protons and neutrons) are described by an effective Dirac equation 
(7 M fc* M — m*)^>(x) = 0, whereas the mesons, which generate the classical mean 
field, are characterized by corresponding covariant equations of motion in the 
Local Density Approximation, with coupling constants that are not density 
dependent. The presence of the hadronic medium modifies the masses and 
momenta of the hadrons, i.e. m* = M + E s (effective masses), k*^ = — E M 
(kinetic momenta). For asymmetric matter the self-energies are different for 
protons and neutrons. This depends on the way the iso-vector mean field is 
described. As discussed in the introduction, in contrast to the iso-scalar case, 
since here one does not have the stringent saturation condition of balancing 
attractive (scalar) and repulsive (vector) contributions, one can apply either 
only the p meson, or both, the p and 5 mesons. We will call the corresponding 
models as NLp and NLpS models, respectively. Furthermore, we will also use 
the non-linear Walecka model without accounting for isospin dependent mean 
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field (TVL-model). 



The parameters used here are summarized in Table 1. The corresponding 
saturation properties of symmetric nuclear matter are E /A(MeV) = —16.0, 
Poifm- 3 ) = 0.16, K(MeV) = 240, m*/M = 0.75 and the symmetry en- 
ergy parameter (for NLp, NLpS) is a 4 (MeV r ) = 30.7 [4,2]. These Lagrangians 
have been already used for flow [25], pion production [5] and isospin tracer 
[13] calculations in relativistic HIC with an overall good agreement to data. 
Moreover it has been recently shown that at high baryon densities the EoS of 
more microscopic DBHF calculations can be well reproduced [5,26]. 

In particular, for the more general NLp5 case one obtains for the self-energies 
of protons and neutrons: 



E a (p, n) = -f a a(p s ) ± fsp s3 



(1) 
(2) 



(upper signs for neutrons), where p s = p sp + p sn , j a = j£ + j",p s3 = ps P - 
Psn, 3% = jp — Jn are the total and isospin scalar densities and currents and 
fa,Lu,p,6 are the coupling constants of the various mesonic fields, see Table 1. 
a(ps) is the solution of the non linear equation for the a field [4,2]. 

Within the relativistic description we expect to see a series of isospin effects 
on the reaction dynamics just due to the changes in the covariant Lorentz 
structure of the isovector interaction. These different isospin contributions 
of relativistic origin can be already seen in the expression for the symmetry 
energy [4,2]: 



J - J sym - — . 



Ikl + 1 



6Et 



m 



fp fs ( ^* 



PB 



(3) 



with Ep = ykp + m* 2 . When including the scalar iso- vector 5-meson the 
asymmetry interaction parameter is given by the combination [f p — fs(p^r) 2 ] of 
the repulsive vector (p) and attractive (5) iso-vector couplings. Thus, in order 
to reproduce the same bulk asymmetry parameter a 4 , we have to increase the 
p-meson coupling when the 5-meson is considered in the iso-vector part of 
the equation of state. The net effect is a stiffer symmetry energy at higher 
baryon densities, due to the j^- quenching of the attractive part, as shown in 

F 

Fig. 1. This is a general feature that has been seen to be important in the 
dynamical simulations [5]. Other relativistic effects of the Lorentz structure 
of the isovector interaction, particularly important for the inelastic channels, 
will be: 



6 



150 - 



100 - 




Fig. 1. Density dependence of the symmetry energy Eq.(3) including only the 
p-meson (dashed line, NL p-model) and both, the p and 5 mesons (dotted line, 
NLp5- model) in the description of asymmetric matter. 

- The vector field S' 4 , Eq.(2), will show a different isospin dependence since 
in the NLpS model the p-meson coupling is larger with respect to that of 
the NLp model. In order to conserve the total energy in inelastic collisions 
we have to include the isospin dependence of the vector fields between the 
ingoing and outgoing channels, as it will be discussed in detail in the next 
section. Thus, the particle production/absorption will be directly sensitive 
to the covariant form of the isovector interaction. 

- Another 5-effect, of interest for transport properties, will be the splitting 
of the neutron/proton effective masses [4,25], see Eq.(4), which will be also 
discussed later on in detail in terms of the energy for subthreshold particle 
production: 

m* = M + Z s (i)=M-f a a(p s )±f sPs3 (+ n, - p) . (4) 



2.2 Hadronic matter in a box 



In order to study hadronic matter, we confine the particles in a cubic box 
with periodic boundary conditions. The number of nucleons and the box di- 
mensions are given as inputs of the simulation and can thus be changed to 
explore different values for the baryon density p# and the asymmetry param- 
eter a of the system. We have fixed the total number of nucleons to ^4=100, 
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neutrons (cc=0) 
neutrons (NLp) 
protons (NLp) 
neutrons (NLpS) 
protons (NLp5) 




E kin (GeV) 

Fig. 2. Kinetic energy spectra {E^ n = E* —m*, p = h\i*) of nucleons for symmetric 
and asymmetric nuclear matter at t = fm/c and fixed baryon density ps = 2.5 po- 
stars: neutron spectra for symmetric matter (a = 0); filled circles and triangles: 
neutron and proton kinetic energy spectra for asymmetric matter (a = 0.2) in the 
NLp model; open circles and triangles: neutron and proton kinetic energy spectra for 
asymmetric matter (a = 0.2) in the NLp5 model. The corresponding solid lines are 
least-square fits according to a Fermi-Dirac distribution at the input temperature 
of T = 60 MeV. 

varying the relative proton and neutron numbers Z and N in order to get 
ol = ttt§ = pBn \ pBp in the range from 0.0 (symmetric nuclear matter) to 0.4 

N+Z PBn+PBp a \ J J 

(which corresponds to extremely asymmetric matter, with N/Z = 2.3). The 
box dimensions are set in order to achieve a baryon density (kept constant 
during the temporal evolution of the system) of roughly 2.5 times the satu- 
ration density, while for the temperature the value of T=60 MeV has been 
chosen. These parameters should resemble the conditions achieved during the 
dense phase in heavy ion collisions at few AGeV energies [27]. 

At t=0 fm/c neutrons and protons are randomly distributed in the box volume 
and their kinetic 3-momenta k* are generated according to relativistic Fermi- 
Dirac distribution functions at the given temperature T: 



1 



l + exp[(£*(k*)-/i*)/T] 



(5) 



In eq. (5) rii (i=n,p) are the occupation numbers for neutrons and protons, 
E* = yk* 2 + m* 2 , while /i* (i=p,n) are the effective chemical potentials for 
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protons and neutrons respectively, given by (e is the energy density): 
de 

— ■ = ^ = p* + Upb T f P PB3 (-n, + p). (6) 



At zero temperature the p* reduce to E* Fi (i=p,n). 

Once both T and a (and hence psi) are fixed, pi, the scalar densities p S i and 
effective masses m* can be calculated self-consistently solving the system of 
equations 



PB^j^yJd^n^k*) (7) 

m* = M — fa(j(p s ) ± f s (ps P - psn) (+ n, - p) (8) 

777* 

d 3 k*->»(k*). (9) 



(2vr)3 J E\ 



For symmetric nuclear matter (a=0) pBi, p S i, Pi an d m* are the same for 
both neutrons and protons so that there is no difference in their momentum 
distributions. In the case of asymmetric matter the momentum distributions 
of neutrons and protons are slightly different due to the different radii of the 
respective Fermi spheres so that p* p ^ p* n ] furthermore, the inclusion of isovec- 
tor scalar meson 5 implies m* ^ m* n . Anyway the initial effective momentum 
distributions of both protons and neutrons just sligthly differ in NL, NLp 
and NLp5 models. This can be seen in Fig. 2, where the nucleon spectra for 
a=0 and a=0.2 are shown. 

We then let the particles (nucleons, A resonances and pions) propagate accord- 
ing to a free kinetic equation plus the collision term, in the spirit of a cascade 
calculation. The mean-field potentials, constant once (ps,o;) are fixed, are ac- 
counted for in the energy-momentum conservation for the various channels 
of the collision term. The boundary conditions are modeled in such a way 
that particles escaping from the box are reinserted at the opposite side with 
the same momentum. In their motion nucleons undergo elastic and inelastic 
binary collisions, thus exciting resonances that further decay into pions, and 
then producing strange particles such as hyperons and kaons. Here we consider 
the lowest mass resonances only, namely A(1232). 

Because of inelastic collisions, a small fraction of the initial kinetic energy of 
the system is lost in particle production and after 150 fm/c nucleons momenta 
are distributed according to a relativistic Fermi-Dirac function with temper- 
ature T=56 MeV. Pion spectra do not appear thermalized yet, though pion 
multiplicities saturate already after t ~ 30 fm/c. Indeed longer time scale are 
needed to reach full thermal equilibrium [28]. 
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For elastic nucleon-nucleon cross sections we use the free parametrizations of 
Cugnon et al. [29], in order to reduce the computing, [30]. For the inelastic 
channels (resonance production) we follow the cross sections of Ref. [31]. Kaon 
production is treated perturbatively, so that it does not affect the overall dy- 
namics. It occurs according to two different creation mechanisms, i.e. baryon 
induced processes such as BB^BYK, and pion induced processes such as 
7rB^YK. In both cases, kaons are created together with Y hyperons (A and 
S) to ensure strangeness conservation and are not reabsorbed by nucleons. 
They can anyway undergo elastic rescattering with nucleons, and thus modify 
their momentum distribution. However, since we are just interested in kaon 
yields and their possible dependence on the symmetry energy, we neglect this 
effect here and do not propagate kaons after production. The elementary cross 
sections for kaon production from BB and nB channels are taken from Refs. 
[32] and [33], respectively. To treat all these scattering processes in infinite 
hadronic asymmetric matter we have modified the original version of the col- 
lision term of the Tubingen RQMD simulation code [18], taking explicitly 
into account the isospin dependence of scalar and vector self-energies S s , 
This turns into an isospin dependence of the threshold of all isospin channels 
and thus affects particle yields, in a way that will be discussed in detail in the 
following section. 



3 Isospin dependent threshold conditions 

The crucial requirement of two-body collisions is the energy-momentum con- 
servation, i.e. the conservation of the total invariant collision energy s, which, 
in terms of ingoing (1,2) and outgoing (3,4) canonical momenta (ki 2 , k^ 4 ), 
reads: 

Sin = (K + fcjOCfyi + k 2^) = (&3 + K)(k 3a + k 4fJ .) = S out . (10) 

If the isospin degree of freedom is not accounted for in inelastic collisions, we 
can simply fulfill energy-momentum conservation using kinetic 4-momenta /c* M 
and the corresponding invariant collision energy 

s* n = (K" + krm, + *y = + krm, + *y = s* out . (n) 

The relation between kinetic and canonical momenta is k*^ = kf — £f, (i — 
1,2,3,4) and the effective masses m* = M + S S j enter in the energy conser- 
vation via the 0-component of the kinetic 4-momentum E* = \/m* 2 + k* 2 . 

The condition (11) is just a constraint on the kinetic plus the effective mass 
energy of the colliding particles and it is sufficient if one considers symmetric 
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(N = Z) matter or restricts itself to elastic collisions or just to the cases where 
the mean field (scalar and vector self-energies) does not change between the 
ingoing and the outgoing particles. In those particular cases Eq. (11) implies 
Eq. (10). 

In the general case of asymmetric nuclear matter, however, the introduction of 
the iso- vector effective mesons causes a splitting in the scalar and in the vector 
fields between protons and neutrons, as discussed in the previous sections (see 
Eqs. (1,2)). For the other hadrons the following expressions in terms of the 
nucleon self-energies will be used. For the four isospin states of the A resonance 
one has [34] (i = scalar, vector) 



E i (A-) = E i (n) (12) 
Ei(A ) = ^(n) + ^(p) (13) 

E i (A + ) = ^I! i (n) + ^i: i (p) (14) 
E 4 (A ++ ) = E,(p) , (15) 

in which the weights (1, |, |) or 'isospin coefficients' are just the square of the 
Clebsch- Gordon coefficients for isospin coupling in the A -<=>- ttN processes. 
In the same framework, for hyperons A, E 0,± we have: 



Ej(A)= 2pM±5M) (16) 
^ )= ! pM+MO) (17) 

E t (E + ) = ^E J (p) (18) 
= ^S,(n) , (19) 

where the hyperon fields are further scaled by a factor | [18], since, according 
to the quark model, there are only two light quarks in a hyperon instead 
of the three light quarks in a nucleon. Pions and kaons are assumed to be 
unaffected by the nuclear mean field and then have zero self-energies. This 
assumption is not strictly correct, at least for kaons, which should be sensitive 
to medium effects and thus shift their production thresholds. Anyway, such a 
shift produced by a kaon potential should be almost the same for both K 0,+ 
[17] and thus it should not significantly affect the K + / K° ratio we will discuss 
here. For this reason we do not take such a potential into account in the 
present work. 
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In inelastic collisions, due to isospin exchange, the scalar and vector self- 
energies between the ingoing and outgoing channels may differ, so that Eq. 
(11) does no longer imply Eq. (10). A typical example is the inelastic process 
nn — > pA~ from which a A~ resonance is formed. It is clear that, even in 
the presence of the p meson only, the vector self-energy of the final channel 
changes and in the general case of accounting for both, the p and S mesons, 
both self-energies (scalar and vector) of the outgoing channel differ from those 
of the ingoing channel. 

In order to properly impose energy-momentum conservation one thus needs 
to replace condition (11) with (10). It then follows from Eq. (10), after some 
algebra, that the effective momenta kg 4 of the outgoing particles (in the cm 
local reference frame in which + kJ; = + k 4 = 0, |k* ut | = |kg| = jk^l) are 
linked to the total energy s in before the collision through the expression 

-(Eg + E°) + y/ 3in + (S 3 + E 4 )2 = ^mf + kf ut + Jmf + kf ut .(20) 



It is worthwhile to stress that in the limiting case of no isospin exchange 
(elastic processes) or isospin independent self-energies (no isospin splittings in 
the scalar and vector self-energies) the general condition (20) reduces to the 
usual energy conservation relation, i.e. s — > \fs* with s* given by Eq. (11). 
Eq. (20) implies the following threshold condition for a given inelastic process: 

s in > (m* 3 + Eg + ml + S°) 2 - (S 3 + S 4 ) 2 . (21) 
.. ' 



This condition reduces to s*„ > (mt 



mz 



whenever we have no isospin de- 



pendence, i.e. for symmetric nuclear matter and in general for the NL model 
and for elastic collisions. The multiplicity of particles produced in the inelas- 
tic processes allowed by Eq. (21) depends on the available energy above the 
threshold, i.e. on the difference As = s in — s , which is affected by the isovec- 
tor channel through both the effective masses and the vector self-energies and 
therefore it changes when considering the NL, NLp and NLpS models. It is 
appropriate to discuss such variations of As for the cases of nn, pp and np col- 
lisions. The rhs of Eq. (21) is not or slighlty affected by both isovector scalar 
(E s ) and vector (E°) mesons depending on the isospin state of the outgoing 
channel. In the particular case of a process nn — ► pA~ (pp — >■ nA ++ ) s is 
respectively given by: 



m p + E a (p) + E°(p) + m A + £ S (A~) + E (A") 



S(p) + S(A-) 
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m n + £ s (n) + E°(n) + m A + E S (A++) + £°(A ++ )] : 
E(n) + £(A++) " 



(22) 



from which it is clear that, according to Eqs.(l-2, 12-15), cancellation occurs 
between Ej(p) and £j(A~). Hence it follows that sq is the same in all the three 
discussed models. For the other inelastic nn and pp channels and in the case 
of np inelastic collisions, which proceed through the excitation of the A or 
A + state of the A resonance, this cancellation is only partial, but in any case 
the net contribution on Sq is moderate. The main effects on As = s in — s are 
therefore coming from Sj„. According to Eq. (10) we have (always in the cm 
local reference frame) 



(£ 1 * + £? + £ 2 * + £°) 2 -(£ 1 + £ 2 r 



(23) 



Since we are dealing with infinite nuclear matter at rest, in which the spatial 
currents vanish, the last term of Eq. (23) is almost negligible also in the 
local frame considered here. Furthermore, due to moderate variations of E* 
induced by m*, £° gives the major contribution to Sj n variations in the different 
interaction models. 



For nn, pp and np collisions we have, respectively 



0\2 



s m = (E* + e; + E^ + Ep) 



(24) 



In the NL model E° is the same for neutrons and protons and hence Sj n is the 
same for all NN collisions (for fixed momenta). When adding the isovector 
vector p-meson only (NLp model), E° increases (see Eq. 2) with respect to 
the NL case, while the opposite occurs for E°. This effect is enhanced in 
the case of NLp5 model, as the introduction of the 5 increases the p-meson 
coupling by about a factor 3. Correspondingly, Sj n increases (decreases) for 
nn (pp) channels when turning the isovector channel on, both without and 
with the isovector scalar meson. At variance, in the case of np collisions due 
to cancellation effects in the sum of (n,p) vector self-energies, Si n is exactly 
the same when the p-meson is included and it will show very small variations 
passing from NLp to NLpS (due to the (n,p) effective mass splitting). 

All that turns out into a gradual increasing (decreasing) multiplicity of A~ 
(A ++ ) states when going from NL to NLp and then to NLp5 model and 
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hence affects the relative populations of different isospin states for pions and 
kaons, as we will widely discuss in the Sect. 5. 



4 Isospin dependent equilibrium conditions 

Since calculations are performed within a box with boundary conditions (par- 
ticles cannot escape) and production and reabsorption processes are consis- 
tently taken into account for A's and pions, we observe that, after a given 
time interval t eq , chemical equilibrium is reached among the different particle 
populations. Indeed, starting from t eq , A and pion multiplicities do not evolve 
anymore, indicating that the production rate equals the reabsorption rate. 
When chemical equilibrium is established, very simple relations hold and the 
chemical potential of produced particles, as pions, can be expressed in terms 
of neutron and proton chemical potentials, p n and p p . In particular, at equi- 
librium, the chemical potential of ir~~ is equal to the difference p n — p p , while 
the chemical potential of ix + is equal to /i p — p, n . Hence one can already make 
analytical predictions for the n~ /n + ratio in the different models used here. 
This ratio can be expressed as: 



where T is the temperature of the system. Neutron and proton chemical po- 
tentials have different values in the models we consider (NL, NLp or NLp5) 
because the self-energies are different. In particular, going from NL to the NLp 
model, we expect that the 7i^/n + ratio increases by a factor exp (4/ p (p Bn — 
PBp)/T), according to the variation of the difference p n — p p (see Eq.(6)). A 
similar, though weaker, increase of the 7i~/n + ratio is expected also when 
going from NLp to NLpS model. 

All these results can be easily fully understood remembering the general rela- 
tion between (n,p) chemical potentials and symmetry energy [2]: 



For a fixed asymmetry a and baryon density ps, one can directly extract the 
equilibrium n~ /n + ratio, Eq.(25). In the NL model (no isovector interaction) 
we have only the weak Fermi contribution to the symmetry energy, that largely 
increases, instead, when the p-coupling is introduced. A further increase, but 
of smaller amplitude, is also observed passing from NLp to NLp5 models, see 



7T /tt+ oc exp (2(/i n - p P )/T), 



(25) 



Pn - Pp = 4:E sym (p B )a. 



(26) 



Fig. 1. 
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Fig. 3. Time dependence of particle multiplicities (as indicated) for symmetric 
(a = 0) hadronic matter at the same fixed conditions as in Fig. 2. The inserted 
panel shows the temporal evolution of kaons in a linear scale. 

Of course these equilibrium considerations do not hold for perturbative kaon 
production, for which reabsorption processes are neglected. 



5 Results 

Here we present the results for excited charge asymmetric matter within the 
effective field models of Sect. 2. The initial conditions of the matter in the 
box have been fixed to a baryon density of ps = 2.5po and a temperature of 
T = 60 MeV, see again Fig. 2. We start from the simpler case of symmetric 
hadronic matter applying only the NL model. Then we will discuss the isospin 
dependence (by varying the iso- vector part of the mean field) of particle pro- 
duction in the asymmetric case beginning with pions and A resonances. A 
similar discussion on strangeness production (i^ + '°) concludes the section. 

Fig. 3 shows the temporal evolution of nucleon, A resonance, pion and kaon 
multiplicities for symmetric matter. Due to isospin symmetry, different isospin 
states for A, 7r and K are almost equally populated. We see that the number of 
nucleons is almost constant throughout the time interval considered and just a 
tiny fraction of them (10~ 2 ) goes into resonances since the first NN collisions. 
Some of these resonances further decay into one-pion channels with typical 
lifetimes of ~1 fm/c and are subsequently reformed in 7rN collisions, according 
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to isospin-dependent cross-sections taken from Ref. [29]. Therefore, the total 
numbers of A and pions grow up in a few fm/c and reach equilibrium values 
of ~3.3 and ~0.4 respectively within 30 fm/c, when formation and absorption 
mechanisms occur with almost the same rate. 

Kaons (K°' + ), on the other hand, do not equilibrate, because they cannot be 
reabsorbed after formation and their number rises almost linearly with time, 
see the inserted panel in Fig. 3. Their multiplicity is indeed well fitted by a 
linear function of time of the form 

N K (t) ~ 1.9 x 1(T 5 t-c 

(c is a constant of the order of 10~ 4 ) where the extremely small rate is due to 
the fact that in the conditions described we are far from the kaon production 
threshold (1.56 GeV for NN collisions in free space) and the observed sub- 
threshold production thus originate from extreme regions of phase space, i.e. 
from the tails of the momentum distributions of the particles involved. Actu- 
ally, kaons could undergo inelastic collisions with other strange particles such 
as hyperons and thus be converted into nonstrange particles. Anyway, studies 
performed at higher energy densities [28] have shown that, even taking such 
processes into account, strange particles require a long time for equilibration, 
of the order of several hundreds fm/c, and hence we can neglect here kaon 
reabsorption as a reasonable approximation, in particular for K°' + . 

The above picture obviously does no longer hold for neutron-rich matter, where 
A _, °, 7r~ and K° are favoured with respect to A +,++ , n + and K + . We will see 
in the following how the yields, and in particular the ratios n~ /ir + and K + /K°, 
will be even sensitive to the Lorentz structure of the isovector interaction. 

5. 1 Isospin effects on pions and A resonances 

Fig. 4 shows the effect of the isovector channels on A and tt multiplicities, 
for the same total density and temperature starting conditions but with a 
fixed asymmetry a = 0.2 (N/Z = 1.5). In the simple NL model the relative 
yields of positively and negatively charged pions and/or the corresponding 
resonances, just follow from the different numbers of nn and pp pairs and, 
consequently, of their inelastic collisions: nn — > pA~ , A~ — > mr~ and pp — > 
nA ++ , A ++ — > pn + . Clearly the ir~ multiplicity increases with asymmetry 
parameter, i.e. as the matter becomes more and more neutron-rich while the 
number of n + accordingly decreases with respect to the symmetric case where 
we have 7r~ = 7r + yields. The further inclusion of a p meson (NLp model), 
stresses such opposite trend in the charged particle production, with the result 
that the former isospin states become even more populated in the case of NLp 
model, while the latter are less populated (solid to dashed lines in Fig. 4). 
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Fig. 4. Time dependence of A resonances (A~° = A" + A ), (A+'++ = A+ + A++) 
and pions (inserted panel), as indicated, for asymmetric (a = 0.2) hadronic matter 
at the same fixed conditions of density and temperature of Fig. 2. Box calculations 
using different models for the mean field are shown: (NL, solid lines) including only 
the iso-scalar part (a, u> mesons), (NLp, dashed lines) accounting for the iso- vector, 
vector p field and (NLp5, dotted lines) taking into account both iso-vector p and 5 
fields. 

Therefore, the number of ir~s, which originate from A~'° decay, also increases, 
while 7r + s decrease (inserted panel in Fig. 4). 



Finally on the basis of the argument discussed in the previous section, one 
expects a further enhancement of the effect by the inclusion of the isovector 
scalar S meson, as confirmed by the results presented in the Figure. 

In Fig. 5 we show the asymmetry dependence of the charged pion multiplic- 
ities. When chemical equilibrium is reached, the 7r~ /n + ratio depends just 
on the difference between neutron and proton chemical potentials and on the 
temperature T. Indeed, if we look at this ratio as a function of the a param- 
eter (right panel) we clearly see the behaviour already noticed from Fig. 4: 
the 7t~/tt + ratio increases when passing from NL to the NLp and then to 
the NLpS model. As explained previously, the n~ /n + ratio should increase by 
a factor exp (4/ p (p j e n — pg p )/T) when passing from NL to NLp model. For 
instance, at the asymmetry a = 0.2, this factor equals 3.6, that approximately 
corresponds to the enhancement observed in the Figure. 

At the same asymmetry, we get for the effective neutron and proton chemical 
potentials: /i* — pf = 34 MeV for NL and NLp models and -16 MeV 
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Fig. 5. Charged pion multiplicities (left) and the ir~ /ir + ratio (right) as a function 
of the asymmetry parameter a at the same conditions of Fig. 2. Statistical errors 
are inside the data points. 

for the NLpS model, due to the nucleoli Dirac mass splitting Eq.(8). Hence, 
considering the change in the difference of effective chemical potentials and 
the value of the constant f p (see Tab. I), we see that the difference p n — p p , 
to which the tt~/tt + ratio is related, changes when going from NLp to NLpS 
model and we expect an increase of this ratio by a factor 1.88. This is also 
nicely confirmed by our results. 



5.1.1 Isospin effects on pion reabsorption 

Since in nuclear collisions equilibrium conditions are likely not reached, it 
is interesting to calculate just the probability of emitting pions in the other 
extreme case, i.e. excluding the pion reabsorption processes. Now the observed 
7r~/7r + ratio should reflect directly the ratio between the cross sections for 7r~ 
and 7r + production. 

As shown in Fig. 6, the tt~ /tt + ratio rises both in the NLp (with respect to the 
NL) model and in the NLpS (with respect to the NLp) model, for any value 
of the asymmetry parameter. In fact, as discussed in the previous sections, the 
addition of each isovector meson (p and S) favours ir~ over tt + production, as a 
result of the higher (lower) energy available in the "parent" nn (pp) collisions. 

We notice that the difference of tt^/ti + ratio observed among the different 
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Fig. 6. Asymmetry dependence of the it /ir + ratio for the same models as in Fig. 
2, but without the inclusion of reabsorption. 

models is less pronounced with respect to the results observed when equilib- 
rium is reached (see Fig. 5). In particular, the results obtained with NLp or 
NLpS are quite close. Moreover, in a nucleus-nucleus collision, pions interact 
with the surrounding matter all along the reaction path. Hence the observed 
final pion multiplicity is not clearly connected to the very first stage of the 
collision and to the high density behaviour of the symmetry energy. In the 
next section we will discuss kaon production, for which we expect a larger 
sensitivity to the models used, since in the conditions considered here we are 
close to the production threshold. Moreover, kaons only loosely interact with 
the nuclear medium and consequently they should carry a better signal of the 
symmetry energy behaviour at large density. 

5.2 Isospin effects on kaon production 

In the case of kaon production the isospin states to be considered are K° 
and K + , which can be produced through a large variety of channels with two 
or three particles in the final state: BB — > BYK, Bit — > YK (B,Y,K 
stand for non-strange baryon, hyperon and kaon, respectively). While, for the 
baryon channels, a proton-proton collision can produce 7r + but not ir~ and a 
neutron-neutron collision, conversely, only 7r~ but not tt + , in the case of kaons 
isospin conservation allows K° and K + to be formed from both initial pp 
and nn scattering processes. However, due to isospin coefficients, K° mainly 
come from nn inelastic channels and K + from pp channels. We thus expect 
their relative yields to change with the asymmetry parameter a, similarly 
to pions without reabsorption, but with a larger sensitivity since kaons are 
produced close to the threshold. Moreover, from the previous discussion on 
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Table 2 

Relative weight of various K°' + production channels for symmetric (a = 0) and 
asymmetric (a = 0.2) matter. The relative width for the AA channels is given in 
units of 10 -4 . Calculations with different symmetry energies are shown too. 
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pion production, also the pionic channels will contribute to isospin effects on 
kaon production. 

The BB — > BYK, Bit — > YK channels can be schematically divided as 
follows: 



(i) NN — ► BYK 

(ii) NA — ► BYK 
(hi) AA — ► BYK 

(iv) irN — >YK 

(v) vrA — ► YK 



where N stands for a nucleon and B for a non-strange baryon. 

At low energy densities, such as those achieved in the present calculations 
and/or in heavy ion collisions at energies below 1.5 AGeV, the pionic channels 
7iN, (iv), give the main contribution [35], together with the NA ones, (ii). 7rA 
and NN collisions contribute just for a few percent to the total K + and/or 
K° yield, while the AA channel is negligible. 

The advantage of a cascade calculation in a box is that we can cleanly evaluate 
the relative weight of the various kaon production channels in the interacting 
matter, i.e. vs. the different choices of the hadron self-energies. Our results are 
reported in Table 2. We see that for all the iso-vector models discussed here we 
are in agreement with previous estimations [35]. We note also, in the case of 
asymmetric matter, the different behavior of the important nN channel when 
we add the isovector mesons: the weight is decreasing for the K + production 
while it is increasing for the K°. This is very promising for our search for fine 
relativistic effects of the isovector interaction on kaon production. 
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Table 3 

Production rate of K +, ° (in units of 10~ 5 ) for asymmetric (a = 0.2) hadronic 
matter for the NL, NLp, NLp5 models. 





dN f t + ^ (xlO- 5 ) 


dN ^P (xlO" 5 ) 


NL 


1.66 


2.18 


NLp 


1.17 


2.91 


NLp5 


1.06 


3.92 



0.006 




time (fm/c) 

Fig. 7. Temporal evolution of K +, ° yields and their ratio (in the inserted panel) for 
asymmetric (a = 0.2) hadronic matter for the same iso-vector models as in Fig. 2. 

Fig. 7 shows the temporal evolution of kaon multiplicities for asymmetric (a = 

0. 2 or N/Z=1.5 corresponding to the asymmetry of 197 Au) hadronic matter 
and their ratio for the NL, NLp and NLpS models. First of all, we observe an 
almost linear time dependence of the kaon yields as already discussed (see Fig. 
2). The isospin effect on the kaon yields essentially turns out in a change of the 
slope, i.e. of the production rate, see Table 3. In particular, for K° (K + ) the 
formation rate increases (decreases) as the symmetry energy becomes stiffer, 

1. e. passing from the NLp to the NLp5 case. 

The whole picture is summarized in Fig. 8, where the kaon yields and their 
ratio are displayed as a function of the asymmetry parameter a. The influence 
of the different iso-vector models obviously increases with a for the K° yield, 
but less for the K°/K + ratio. In fact, as already shown in Fig. 7 and Table 3 
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Fig. 8. Asymmetry dependence of the K +, ° yields (left) and of the K + /K° ratio 
(right) for the different iso-vector models as in Fig. 2. Statistical error bars are also 
shown, when larger than the data points. 

for a fixed asymmetry, the K + yield, which decreases of about 25% and 10% in 
NLp and NLpS models respectively, presents an isovector meson effect slightly 
weaker than the K° yield, which exhibits increasing variations of about 35% 
for both NLp and NLpS models. 



Generally the isospin effects on kaon yields originate from at least two different 
mechanisms: the isospin dependencies of the yields and of the threshold 
energies. We will now discuss in detail these two effects. 

Since kaon production from pionic channels (which give the main contribution) 
is competitive with pion reabsorption, it occurs before secondary ttN — > A 
processes set in. Therefore, K° and K + kaon yields are mainly affected by the 
7r~ and 7r + abundances before reabsorption, see Fig. 6. In particular, as K° 
(K + ) mainly come from tc~N (tt + N) processes, the enhancement (reduction) 
of the 7r~ (7r + ) multiplicities in the NLp5 model is partially responsible for a 
corresponding increase (decrease) of the K° (K + ) yields. 

The isospin effect on the kaon production threshold is more complicated, due 
to the many isospin channels that may affect the threshold (and hence the 
yields) differently. However, contributions in opposite directions originating 
from different channels do not exactly cancel each other. In order to clarify 
the observed influences of the models on the yields in a transparent way, it 
is appropriate to focus on some demonstrative channels. Such examples are 
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the isospin channels rm — > E K° and pn + — > T, + K + for K°' + production, 
respectively. 

ttN channels offer the advantage of two particles in the final state (instead of 
the three that we have in the BB case); moreover, since also kaons (and not 
only pions) are, in the present calculations, unaffected by nuclear mean field, 
Eqs. (20-23) are simplified (E i2 = E i7r = Ej 4 = E^ = 0). 

According to Eq. (20) the corresponding threshold conditions read 

Sin{nn-) > (M s - + M K o + E sS - + E E _) 2 + S|_ (27) 
Sin(pn + ) > (M s+ +M K+ + E sS+ + E s+ ) 2 + E| + , (28) 

where M s - + M K o = 1695Ge\/ and M s+ + M K+ = l6S3GeV. 

As discussed in the previous section, for collisions involving neutrons (protons) 
Si n increases (decreases) when going from NL to NLp and then to NLpS 
model. Therefore, the available energy for K° production increases as the 
symmetry energy becomes stiffer, leading to an enhanced K° yield. Obviously 
the trend for the isospin behavior of the K + yield will be opposite. 



6 Conclusions 

We have investigated the high density behavior of the iso-vector of the nu- 
clear interaction, which is still poorly known experimentally, controversially 
predicted by theory, but of great interest in extreme nuclear systems. Here we 
discuss its dynamical effects on particle production for an idealized system, the 
asymmetric hadronic matter. In this respect relativistic cascade calculations 
have been performed in a box with periodic boundary conditions. The advan- 
tage is that in this way we can cleanly follow all the contributions of the various 
isospin channels to the production and reabsorption, and the corresponding 
dependence on the effective field structure of the isovector interaction. 

The paper represents a detailed study, in a relativistic frame, of the produc- 
tion/absorption mechanisms of energetic particles (pions and kaons) in excited 
nuclear matter at high isospin and baryon density and temperature. These 
are transient conditions in realistic charge asymmetric heavy ion collisions at 
intermediate energies. The yields are largely influenced by the isospin depen- 
dence of the nucleon self-energies. In particular, a given isospin state of the 
produced particles is differently populated varying the high density behaviour 
of the symmetry energy. 

It has been found that, after ps 30 fm/c, chemical equilibrium is reached. The 
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7t~/tt + ratio depends on the model considered, according to the difference 
between neutron and proton chemical potentials. Excluding the reabsorption 
process, in order to measure just the pion emission probability in situations 
when equilibrium is not necessarily reached, we find a less pronounced de- 
pendence of the 7r~/7r + ratio on the symmetry energy behaviour. This seems 
somehow confirmed by theoretical transport calculations on Au+Au collisions 
[5,34]. 

The isovector meson effect on kaon rates appears larger than on pions (exclud- 
ing reabsorption) because kaons are produced close to the threshold. It should 
be noticed that low energy kaons do not suffer from secondary reabsorption 
effects neither in infinite hadronic matter nor in collisions of finite nuclei due 
to strangeness conservation. Hence, in a nuclear collision, they are expected 
to better signal the behaviour of the symmetry energy during the first stage 
of the collision, when they are formed. As an important result, the K° and 
K + yields are oppositely affected leading to a crucial isospin dependence of 
the K + /K° ratio, which one would expect to observe in heavy ion collisions 
at subthreshold energies. Such an investigation is being performed and will be 
discussed in a future work. 

In conclusion the results presented here appear very promising for the pos- 
sibility of extracting information from terrestrial laboratories on the Lorentz 
structure of the isovector nuclear interaction in a medium at densities of astro- 
physical interest. The use of radioactive beams at relativistic energies would 
be extremely important. 
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